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THERMAL  RESPONSE  OF  GRAPHITE  EPOXY  COMPOSITE 
SUBJECTED  TO  RAPID  KEATING 


1.  INTRODUCTION 

Due  to  their  high  strength  to  weight  ratio  at  low  temperatures,  fiber  reinforced  organic  matrix 
composites  have  been  successfully  utilized  in  numerous  military  and  space  applications.  However, 
intense  heating  produced  by  fire  or  laser  or  irradiation  can  adversely  effect  the  integrity  of  composite 
structures  due  to  the  rapid  degradation  of  mechanical  properties  (tensile  and  shear  strength,  elastic 
moduli,  etc.)  at  temperatures  above  200°C  (1).  Severe  thermal  irradiation  may  also  induce  undesirable 
geometric  changes  through  ablation  of  critical  load  carrying  members  or  by  generation  of  high  stress 
concentration  due  to  localized  burnthrough  (2).  Thus,  accurate  modeling  of  thermal  response  consti¬ 
tutes  an  important  element  in  assessing  the  overall  reliability  of  composite  structures  subjected  to  rapid, 
high-intensity  heating. 


In  the  present  study,  a  simplified  one-dimensional  mathematical  approach  is  presented  for  tran¬ 
sient  heat  transfer  analysis  of  composite  plates  exposed  to  surface  heating  This  model  is  then 
employed,  in  conjunction  with  an  experimental  stydy,  to  verify  recently  proposed  (3)  thermophysical 
properties  of  AS/3501-6  graphite  epoxy  laminate.  In  a  separate  investigation,  computational  results 
(temperature  distributions  and  ablative  characteristics)  from  the  thermal  analysis  will,  be  coupled  with 
detailed  stress  analyses  to  evolve  a  thermomechanical  failure  criterion  for  AS/3501-6  laminates  under 


combined  laser  irradiation  and  uniaxial  tensile  loading. 

*  • 

2.  ANALYTICAL  DEVELOPMENT 

2.1  Heat  Transfer  Equations  for  One-Dimensional  Model 

Since  the  current  study  is  concerned  largely  with  scenerios  where  a  uniform  high-intensity  energy 

' 

flux  is  applied  over  an  area  whose  characteristic  dimensions  are  significantly  greater  than  the  composite 
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thickness,  a  one-dimensional  approach  is  adopted.  Keat  transfer  is  assumed  to  occur  by  conduction  in 
the  thickness  direction,  i.e.,  aormal  to  the  plane  of  the  fibers.  Thermal  decomposition  of  the  organic 
matrix  material  (e.g.,  epoxy  resin)  is  approximated  (3)  by  elevation  of  the  bulk  heat  capacity  over  an 
appropriate  temperature  range;  thus,  temperature  dependent  material  properties  are  included  in  the 
treatment.  In  addition,  the  fiber  reinforcing  component  (e.g.,  graphite  or  boron)  as  well  as  any  residual 
matrix  constituents  are  assumed  to  undergo  a  simple  sublimation  reaction  when  the  surface  temperature 
reaches  a  predetermined  constant  value.  Mathematically  this  surface  ablation  may  be  considered 
equivalent  to  a  classical  melting  problem  (4)  in  which  the  melt  is  instantaneously  removed.  For  subse¬ 
quent  numerical  convenience,  a  convective  coordinate  system  (5.6)  is  adopted  which  moves  continu¬ 
ously  with  the  ablating  fibers  as  indicated  in  Fig.l,  i.e., 

z  -  X  -  xb,  (1) 

where  z  and  x  are  the  convected  and  global  coordinates,  respectively,  and  xh(t)  denotes  the  instantane¬ 
ous  position  of  the  receding  surface.  With  respect  to  the  moving  coordinate,  the  one-dimensional 
Fourier  conduction  equation  is  given  by 

+  „C,K— -PC,— ,  (2) 

where  T(,z,t),  k.p  and  Cp  represent  the  temperature,  thermal  conductivity,  density,  and  heat  capacity, 
respectively.  The  quantity  V  =  dxjdt  denotes  the  surface  recession  rate  and  is  regarded  as  an  unk¬ 
nown  function  of  time. 

Subsequent  to  the  onset  of  abiation  (r  <  t, >),  boundary  conditions  at  the  irradiated  surface 
(z  -  0)  include: 

no,/)  —  T,  -  constant  (3) 

and  a  straightforward  energy  balance  expressed  by 

+  /,  +  /„  (4) 

in  which  H,  is  an  effective  heat  of  ablation  (assumed  constant)  and  /0  represents  the  incident  heat  flux. 
Surface  energy  losses  due  to  radiation  and  convection  are  I,  and  le,  respectively,  in  Eq.  (4),  and  the 
coupling  coefficient  a  denotes  the  fraction  of  the  incident  energy  which  is  absorbed  by  the  receiver. 
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The  radiation  toss  cited  above  is  computed  according  to 

/,  -  <rt  (T$  -  T*).  (5) 

where  cr  is  the  Stefan-Boltzmann  constant;  <  is  the  surface  emissivity;  and  Ta  represents  the  surround¬ 
ing  air  temperature. 

The  convective  heat  loss,  Ic,  resulting  from  the  flow  of  high  velocity  air  parallel  to  the  irradiated 
surface  is  given  by  the  Newton  heat  transfer  equation,  viz., 

Ic-h(T,-T,).  (6) 

where  h  and  Tr  are  the  convection  coefficient  and  air  recovery  temperature,  respectively.  In  the  present 

study  a  laminar  boundary  layer  is  considered  in  which  case  the  Pohlhaussen  flat  plate  equation  may  be 

used  to  estimate  the  convection  coefficient.  A  simplified  form  (7)  of  this  equation  is: 

h  (btu/ia1  -  sec  -  T)  -  0.3—  Re^l  (7) 

5 

where  ka  is  the  thermal  conductivity  of  air,  s  represents  the  distance  from  the  leading  edge  of  the  boun¬ 
dary  layer,  and  Re  is  the  Reynolds  Number.* 

Due  to  the  low  transverse  conductivity  exhibited  by  many  organic  matrix  composites  at  elevated 
temperatures,  a  relatively  small  temperature  rise  is  frequently  observed  at  the  unirradiated  surface  over 
the  exposure  times  of  interest.  Consequently  the  rear  surface  is  assumed  to  be  insulated,  i.e., 

~  <■  0  at  z  «•  I  —  (8) 

where  /  is  the  initial  laminate  thickness. 

It  may  be  noted  that  prior  to  ablation  (r  <  f0.  T(0.t)  <  T,),  Eqs.  (1)  and  (2)  are  appi:cable  with 
xb  —  V  *■  0.  Appropriate  boundary  conditions  are  then  expressed  by  Eqs.(3-8)  with  V  »  0  and  the  sur¬ 
face  temperature  T,  regarded  as  an  unknown  function  of  time. 


*  The  semi-empirical  expressions  *tven  in  Ref.  (7)  were  utilized  io  evaluate  the  dependence  of  It,  and  Re  on  surface  tempera¬ 
ture  xiso.  the  thermodynamic  analysis  in  Ref.  17)  was  employed  to  compute  T,  <Eq.  (6))  in  terms  of  the  airstream  velocity  and 
Tr,.  this  approach  assumes  a  constant  value  of  the  recovery  factor  equal  to  0.88. 
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2.2  Numerical  Solution 

Solutions  to  the  nonlinear  boundary  value  problem  defined  by  Eqs.  (1-8)  were  obtained  numeri¬ 
cally  using  a  modified  Crank-Nicholson  finite  difference  method.  Figure  l  gives  a  schematic  representa¬ 
tion  of  the  one-dimensional  array  of  nodal  points  emploved  in  the  calculations.  To  efficiently 
characterize  the  steep  temperature  gradient  near  the  irradiated  surface,  a  variable  mesh  spacing  is  used 
in  which  the  the  nodal  separation  increases  with  increasing  distance  from  the  front  surface.  Since  Eqs. 
(1-8)  are  written  with  respect  to  the  convected  coordinate  z,  the  entire  finite  difference  grid  is  uni¬ 
formly  translated  along  with  the  receding  front  surface  as  the  solution  proceeds.  However,  because  the 
global  coordinate  of  the  unirradiated  surface  (x  —  /)  is  fixed,  the  total  number  of  nodes  actually  utilized 
decreases  as  ablation  continues;  i.e.,  once  a  node  has  traversed  the  rear  surface  it  is  ignored  in  subse¬ 
quent  calculations. 


With  respect  to  the  generic  sequence  of  interior  nodal  points  i  -  1.  /.  i  +  1  (rig.  1),  finite 
difference  approximations  for  the  first  and  second  temperature  derivatives  at  node  i  are: 

9T  —  Ti*u  +  DTj.y  ~  . 

dz  tJ  r(r  +  l)(z,  -  z,_i)  8 

and 

d2r]  2(rr,_u  +  Tj+u  ~  (r  +  1)TJ 

_  ■*  *  .  ,  .  ,  .  •»  »  w  O/ 

9z2  ,j  r(r  +  l)(z,  -  z,_t)2 

where  r  s  (r,+i  -  z,)/(z,  -  z,_t)  and  the  subscripts  i  and  J  on  T and  its  derivatives  indicate  evaluation 
at  coordinate  z,  at  time  t,.  As  described  in  tlef.  (8)  the  Crank-Nicholson  method  entails  using  in  Eqs. 
(1-8)  the  following  average  spatial  derivatives  over  a  generic  time  increment  It  a  ty+1  -  i,. 

U  .  [il 

d2  >.m  “  ** 

and 
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where  the  subscript  m  implies  evaluation  at  the  mean  time  over  the  increment  tm  —  (r,+1  +  t,)/2.  The 
corresponding  average  time  derivative  over  the  interval  It  is  represented  by  the  forward  difference 
expression: 

3 T  _  ~  Tt.j 

Substitution  of  Eqs.  (9),  (10),  and  (11)  into  Eq.  (2)  yields  the  general  difference  equation  to  be 
satisfied  at  all  nodal  points,  viz.. 


1  pCf(z,  — ,  z,_i)J 

2  2(z,  -  z,_ |)  ' 

(dk/dT) \t 

3  4r2(l  +  r)2(z,  -  Z,_t)JpC, 

The  material  properties  (p,C,.k.dkldT )  referred  to  in  Eq.  (12)  are  generally  temperature  dependent 
and  are  to  be  evaluated  at  the  mean  temperature  T,im  existing  at  node  i  over  the  time  It.  Similarly,  the 
surface  recession  rate  Vm  appearing  in  X  2  is  an  average  value  associated  with  time  tm. 


According  to  Eq.  (3)  one  boundary  condition  at  the  inadiated  surface  (node  1  in  Fig.  1)  is  simply 
the  constant  temperature  condition 

7r(0,tj)  a  Tu  -  T„  (13) 

which  is  valid  for  t  >  t^  The  conservation  of  energy  condition,  Eq.  (4),  is  also  applicable  at  node  l 

and  may  be  conveniently  expressed  in  difference  form  by  introduction  of  a  virtual  node  at  z  •  ~Zi  as 

indicated  in  Fig.  1.  A  straightforward  substitution  of  Eqs.  (9a)  and  (10a)  into  Eq.  (4)  then  gives 

7Y,*.  +  7Y,-  To.,*,-  (14) 
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in  which  Q  3  a  /0  +  /,  +  /,.  is  constant  by  virtue  of  Eq.  (3)  providing  IQ  does  not  vary  with  time.  By 
applying  the  general  expression,  Eq.  (12),  at  node  1  and  combining  the  result  with  Eq.  (14),  the  vir¬ 


tual  temperatures  (r0y+i  and  r0y)  are  eliminated.  This  gives 

,r-r  pH.V.-Q  \„C,2lV. 

2  *i.j+\  *2j  "  k  |  k - + 

where  all  material  properties  are  evaluated  at  the  surface  temperatue  Ts. 


( dk/dT )  Z2 


The  zero  heat  flux  condition  at  the  unirradiated  surface  is  implemented  numerically  by  again  using 
a  virtual  node  at  z  —  2x,  -  zf_,  (see  Fig.  1).  A  simple  two-term  series  expansion  about  node  r  is 
employed  to  estimate  the  temperature  gradient  at  the  insulated  surface,  viz., 

¥  -¥  *--»•  <“> 

tmtf.m  ®2  r.m  r.m 

where  the  second  equality  follows  from  Eq.  (8),  and  ft*  refers  to  the  temporally  averaged  distance 
between  node  r  and  the  unirradiated  surface  over  the  time  increment  &r.  The  quantity  ft*  can  be 
expressed  by  the  recursive  relation: 

ft  ,+i  +  ft :  Vm  Ar 

»*-  ;  '2-  ;  -ft;- -y-.  (17) 

in  which  8y  is  the  separation  existing  at  the  begining  of  the  time  increment,  i.e.,  at  t  -  t,.  Combination 
of  Eqs.  (9),  (10),  (16),  and  (17)  yields: 


2(7-,,*,  +  TrJ)kt-  (r,.u+|  +  7V_ly)<ft4-  1/2)-  (rf+l./+,  +  T,+1,)(R4  +  l/2).  (18) 

where 

28  y  -  vm\t 

4”  2<zf-  .%_,)’ 

Finally  the  temperatures  T,* ly>i  and  y  at  the  virtual  node  can  be  eliminated  by  invoicing  the  gen¬ 
eral  difference  relationship  (Eq.  (12))  at  node  r  and  combining  the  result  with  Eq.  (18).  This  pro¬ 
cedure  provides: 


-ft  i  +  Xj  + 


(ftl  +ft;)(2ft«-  1) 
2ft«  +  1 


4ft4(ft,+ftj) 

Tr-u+x  +  2(1  +  ft,)  -  — ^  t  :  rr,+l  - 


r  (ft  ,  +  ft  ;)(2ft4  -  1)  4ft4(ft,+ftt), 

P  2vTT1  r'-'-/  +  [2(l  "ft,)+  >4+l  T" 
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Equations  (12),  (IS)  and  (19)  define  an  implicit  numerical  method  for  solution  to  the  present 
boundary  value  problem;  given  the  solution  at  a  generic  time  tJy  simultaneous  solution  of  these  equa¬ 
tions  provides  the  temperature  distribution  and  recession  rate  at  ti+l.  Since  the  difference  equations  are 
obviously  nonlinear,  their  solution  at  each  time  step  was  obtained  iteratively  using  the  method  of  suc¬ 
cessive  approximations.  For  a  typical  time  increment  Ar  -  ry+1  -  f,,  the  following  algorithm  was 
adopted.  In  the  first  iteration,  the  quadratic  terms  appearing  in  Eqs.  (12)  and  (19)  as  well  as  the  values 
of  \|,A2,Aj,  and  \4  were  evaluated  using  the  known  nodal  temperatures  and  recession  rate  existing  at 
tj.  The  resulting  linear  versions  of  Eqs.  (12)  and  (19)  were  solved  simultaneously  for  the  temperatures 
of  nodes  2  through  r  at  ry4l.  Then,  using  Eq.  (15)  in  conjunction  with  the  recently  computed  Tjj+u  30 
initial  estimate  of  Vm  was  calculated.  The  process  was  then  repeated  using  the  most  recent  tempera¬ 
tures  and  Vm  to  generate  a  new  (updated)  linear  form  of  Eqns.  (12)  and  (19).  At  each  time  step,  the 
calculations  were  continued  until  uniform  values  of  Vm  (within  two  percent)  were  obtained  in  succes¬ 
sive  iterations. 

As  suggested  in  Sec.  2.1  finite  difference  exp.essions  analogous  to  Eqs.  (12),  (15),  and  (19)  may 
be  developed  for  pre-ablation  regime  (r  <  to)  by  noting  that  the  surface  temperature  TX  )  is  variable 
and  that  Vm  —  0.  The  resulting  simultaneous  difference  equations  are  again  nonlinear  necessitating  use 
of  an  iterative  procedure.  A  method  similar  to  that  described  above  was  employed  in  which  conver¬ 
gence  was  based  on  constancy  of  surface  temperature  (rather  than  Vm)  in  successive  iterations. 

EXPERIMENTAL  PROCEDURE  AND  MATERIAL  PROPERTIES 

3.1  Laser  Irradiation  Tests 

To  evaluate  the  analytical  model,  laminated  composite  panels  were  subjected  to  rapid  heating 
using  the  15  kw,  continuous  wave,  COj  laser  at  NRL.  Twenty-ply  AS/3501-6  graphite  epoxy  coupons 
having  in-plane  dimensions  5.6  x  10  cm.  were  fabricated  by  the  NRL  Chemistry  Division.  The  panel 
thickness  was  2.54  mm.  and  the  piy  layup  sequence  was 

[  +  45,  -  45.  90.  0.  0.  0.-  45.  +  45.  0.  0.  ],. 
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where  the  subscript  denotes  symmetry  about  the  laminate  midthickness  and  0.  +  45.  —  45,  etc.  indicate 
the  angular  orientation  of  the  lamina  fibers  with  respect  to  the  lengthwise  direction  of  the  specimen. 
To  measure  the  thermal  response  during  irradiation,  one  thermocouple  was  embedded  in  the  11th  ply 
during  specimen  fabrication  and  a  second  was  placed  ?\  the  rear  surface  (ply  20).  In  addition,  an  optical 
pyrometer  was  employed  to  establish  the  front  surface  temperature.  Testing  was  conducted  at  several 
intensity  levels  (/<>)  which  were  generated  using  a  fixed  25mm  beam  diameter  with  variable  power  out¬ 
put.  To  simulate  aerodynamic  cooling  effects,  a  Mach  0.3  airflow  was  applied  parallel  to  the  irradiated 
surface  during  each  test.* 

3.2  Thennophysical  Properties 

The  temperature  dependence  of  the  material  properties  of  AS/3501-6  which  were  employed  in  the 
numerical  calculations  are  indicated  in  Fig.  2.  These  properties  were  obtained  from  Ref.  (3)  with  minor 
modifications.  The  density  variation  p(D  proposed  in  Ref.  (3)  was  increased  uniformly  by  6  percent 
based  on  a  higher  room  temperature  value  observed  for  the  composite  used  in  the  present  experiments. 
The  plateau  in  the  heat  capacity  curve  CP(T)  for  343  <  T  <  510°C represents  additional  heat  absorbed 
during  resin  (matrix)  decomposition;  effective  heat  capacity  values  over  this  range  were  computed  (3) 
using  a  theoretically  derived  latent  heat  of  9%  J  per  gram  of  resin.  In  the  present  work  the  abrupt  vari¬ 
ations  in  Cp  at  343  and  510*C  suggested  in  Ref.  (3)  were  replaced  by  more  gradual  changes  occurring 
over  28°C  intervals.  Similarly,  a  smooth  variation  in  conductivity  k(T)  and  a  less  severe  28°C  transi¬ 
tion  in  in  p(T)  were  introduced  at  538’C.  These  modest  shape  changes  in  the  property  curves  pro¬ 
duced  no  significant  change  in  analytical  results  and  were  adopted  soley  to  provide  enhanced  stability 
and  convergence  characteristics  in  the  numerical  calculations. 

As  noted  in  the  Cf(T)  curve  of  Fig.  2,  the  values  of  Ht  and  Ts  for  graphite  sublimation  were 
assumed  to  be  43  kJ/g  and  3316’C,  respectively.  These  properties  were  also  obtained  from  Ref:  (3). 


•  For  ihese  experiment*  the  appropriate  value  of  s  to  be  used  in  the  computation  of  the  convection  coefficient,  Eq.  <7),  is  10.48 

cm. 
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4.  DISCUSSION  OF  RESULTS 

The  accuracy  of  the  analytical  procedure  and  assumed  material  properties  may  be  assessed  by 
comparison  of  the  experimental  and  numerical  results  shown  is  Figs.  3,  4,  and  5  for  incident  beam 
intensities  (/o)  of  0.232,  1.33,  and  2.79  kw/cm2,  respectively.*  It  is  noted  that  the  calculated  and  meas¬ 
ured  temperatures  agree  reasonably  well,  particularly  at  the  lowest  intensity  level  (Fig.  3)  where  the 
disparity  is  less  than  45°C  over  most  of  the  exposure  time  investigated.  Examination  of  the  numerical 
results  for  the  2.79  kw/cm2  case.  Fig.  5,  reveals  that  the  rapid  increase  in  interior  (ply  ID  temperature 
after  approximately  4  seconds  is  due  to  the  close  proximity  of  the  advancing  irradiated  surface  to  the 
(global)  x-  1.4  mm  coordinate.  The  comparable  increase  in  the  embedded  thermocouple  response 
suggests  that  the  ablation  characteristics  of  the  laminate  are  modeled  quite  well. 

Figures  3-5  also  indicate  good  correspondence  between  predicted  and  measured  front  surface  tem¬ 
peratures.  The  lower  temperature,  1860°C,  measured  after  irradiation  for  9.2  sec  at  0.232  kw/cm2  is 
reasonably  close  to  the  corresponding  analytical  result  of  2132°C.  No  ablation  is  predicted  for  this  test 
condition.  On  the  other  hand,  significant  ablation  is  projected  for  the  1.33  and,  2.79  kw/cm2  beam 
intensities  and  the  measured  surface  temperatures,  viz.  2938  and  3371*C,  respectively,  compare  very 
well  with  the  assumed  3316°C  temperature  level  necessary  to  sustain  graphite  sublimation. 

Figure  6  compares  the  analytical  and  experimental  mass  loss  as  a  function  of  exposure  time  for 
material  subjected  to  an  average  intensity  of  1.81  kw/cm2.  The  analytical  lm(t)  was  determined  from 
the  easily  derived  expression 

±m(t)  -  A  J^(f)  ( po  -  p(.T)]dx, 

whrtre  p(D  is  given  in  Fig.  2;  po  represents  the  initial  (room  temperature)  density,  and  A  is  the  total 
beam  area.  The  temperature  distribution  T  —  T(x,t)  as  well  as  the  instantaneous  position  of  the  irradi¬ 
ated  surface  xt(r)  were  obtained  from  the  numerical  solution.  Although  excellent  agreement  between 


*  These  intensities  refer  to  maximum  values  existing  at  the  center  of  the  beam.  Also,  the  surface  emissivity  (a)  and  absorption 
coefficient  (al  were  both  asngned  a  value  of  9.92  in  the  numerics!  calculations. 
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predicted  and  measured  mass  loss  is  evident,  it  should  be  noted  that  an  average  rather  than  peak  inten¬ 
sity  was  employed  in  the  computation  of  T(x,t)  and  x0.  Use  of  the  peak  flux  (2.57  kw/cm2)  in  the  cal¬ 
culations  would  lead  to  significantly  greater  mass  loss  rates.  However,  utilization  of  the  average  inten¬ 
sity  appears  reasonable  since  the  experiments  themselves  represent  a  spatially  average  response  to  irra¬ 
diation  which  includes  the  effects  of  non-unifo.m  beam  intensity. 

Figures  7  and  8  show  the  computed  variation  in  recession  rate  V  with  time  for  /0  values  of  1.81 
and  2.79  kw/cm2,  respectively.  The  inset  curve  in  each  figure  indicates  the  short-time  response.  Once 
the  sublimation  condition  is  attained  at  t  —  (0,  an  initial  rapid  increase  in  V  is  apparent  followed  by  a 
significant  period  during  which  the  irradiated  surface  advances  at  a  relatively  uniform  velocity.  Figure  7 
also  indicates  that  at  longer  exposure  times  V  again  increases  somewhat  as  the  moving  front  surface 
approaches  the  back  face.  This  latter  effect  simply  reflects  a  less  severe  temperature  gradient  ahead  of 
the  front  surface  which  is  caused  by  the  insulation  condition  imposed  at  x  -  /.  This  reduced  gradient 
implies  that  the  second  term  on  the  right  hand  side  of  Eq.  (4)  becomes  larger  (less  negative)  indicating 
higher  \alues  of  V. 

The  nearly  constant  recession  rate  displayed  over  much  of  the  exposure  time  may  be  approxi¬ 
mated  by  consideration  of  the  steady  state  solution  of  Eqs.  (1-7)  for  a  semi-infinite  plate.  Regarding  V 
and  all  material  properties  as  constant  and  requiring  that  T  —  To  as  z  -»«•’,  the  solution  is 

Hz)  -  (Ts  -  T0)  exp  +  To. 

which  may  be  inserted  into  Eq.  (4)  to  give 

v  a/p  +  l,  +  fc 

m  PIH,  +  C,iT,  -  Toll- 

Evaluation  of  /,,  lc,  and  material  properties  at  Ts  —  3316  TT  provides  the  following  dependence  of  V  on 
Io  for  the  present  experiments: 

K(mm/ sec)  -  0.156  /0  -  0.150. 

where  /0  has  units  of  kw/cm2.  As  indicated  by  the  dashed  lines  in  Figs.  (7)  and  (8) .  this  relationship 
provides  a  good  estimate  of  the  average  velocity  oyer  the  exposure  times  indicated  and  may  prove  use¬ 
ful  in  projecting  plate  bumthrough  times  for  AS/3501-6  composite. 
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Although  the  numerical  treatment  presented  in  this  study  is  adequate  for  most  engineering  studies 
associated  with  elevated  temperature  composite  behavior,  improved  modeling  is  undoubtedly  possible. 
Further  examination  of  the  calculated  thermal  responses  shown  in  Figs.  3-5  indicates  that  the  analysis 
consistently  underestimates  the  measured  interior  temperatures  and  tends  to  overestimate  the  rear  sur¬ 
face  values.  Improved  correlations  at  the  back  surface  may  result  if  reduced  thermal  conductivity 
values  were  introduced  at  lew  tempeteratures,  and  Ref.  (9)  contains  experimental  conductivity  data 
which  would  support  such  an  approach.  In  addition,  a  recent  study  (1)  of  graphite  epoxy  exposed  to 
flame  heating  has  established  that  exothermic  resin  decomposition  is  possible  and  a  kinematic  expression 
quantifying  the  associated  heat  generation  has  been  proposed.  Inclusion  of  such  a  heat  source  term  in 
the  conduction  equation  should  lead  to  higher  predicted  temperatures  inside  the  laminate  and  better 
agreement  with  experimental  results. 

5.  SUMMARY 

The  primary  conclusions  from  the  present  study  of  the  thermal  response  of  graphite  epoxy  com¬ 
posite  subjected  to  rapid  heating  include: 

•  The  proposed  one-dimensional  numerical  model  coupled  wna  the  material  properties  contained  in 
Ref.  (3)  provides  temperature  distributions  and  ablative  response  which  are  in  good  agreement 
with  experimental  results  from  laser  irradiation  tests.  When  incorporated  with  appropriate  stress 
computations,  the  present  analysis  provides  a  practical  method  for  projecting  the  effects  of  severe 
thermal  environments  on  composite  plate  structures. 

•  The  present  numerical  results  indicate  that  a  steady  state  thermal  response  is  approximately 
achieved  over  a  significant  portion  of  the  exposure  times  investigated.  Analytic  expressions  for 
the  temperature  distribution  and  recession  rate  under  these  conditions  have  been  given. 

•  Suggestions  for  improved  modeling  include  incorporation  of  exothermic  matrix  decomposition 
kinetics  and  the  inclusion  of  reduced  bulk  thermal  conductivity  values  at  low  temperature. 
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Figure  1  —  Convected  Coordinate  System  and  Finite  Difference  Grid  Employed 
in  One- Dimensional  Analysis. 
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Figure  3  —  Comperaon  of  Measured  and  Predicted  Temperature  Response  for  AS/3501-0  Uminete. 

Beam  Intensity  «  0.232  lew/ cm-. 


TEMPERATURE  (°C) 


GRIFFIS  MASUMURA.  AND  CHANG 


Figure  4  —  Comparison  of  Measured  and  Predicted  Temperature  Response  for  AS/3501-6  Laminate. 

Beam  Intensity  —  1.33  kw/cm*. 
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TIME  (sec) 

Figure  5  —  Compel  non  of  Measured  and  Predicted  Temperature  Response  for  AS/3501-6  Laminate. 

Beam  Intensity  **  2.79  kw/cm*. 
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Figure  6  —  Experimental  and  Aroiytxaily  Derived  Maae  Loee  u  a  Function  of  Exposure  Time. 
Avera«e  Beem  Intenaity  “  LSI  IcWcnr. 
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RECESSION  RATE,  V(mm/tec) 
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F«ure  I  —  Variation  of  Rmmmm  Ram  with  Exposure  run*  for  2.79  Irw/crtr1  Beam  Intensity. 
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